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ABSTRACT 


Subscripts: 


Ensieo’etder oempling was used to oliminato the 
integral haimonics from tlie flutter spectra corre- 
sponding to a case-mounted static pressure trans- 
ducer. Using tie optical displacemont data, it was 
demonstrated that the blade-order sampling of pressure 
data may yield erroneous results duo to the inter- 
ference caused by blade vibration. Two methods are 
proaenced which effectively oliminato this Incor- 
feronco yielding the blado-preaouro-difforonce spec- 
tra. The phase difforanco between the dif f erential- 
preasure and the displacement spectra was ovnluatodc 


f flutter 

i Corresponds to i'th nodal dlamoter 
L loft 
P pressure 
R right 
r rptational 


NOMENCLATURE 


Superscripts: 


A amplitude 


D displacement 


D diaplaccment in time units or in IE periods, E"1 


P pressure 


d displacement (cm), d “ D*U (Fig. 1) 


INTRODUCTION 


E engine order 

j imaginary unit 

n number of points for spectral analysis 

N nodal diameter 

P unsteady pressure 

P ateady-stato pressure 

r radi.ll coordinate (Fig. 1) 

T sampling period 

t time 

U wheel speed 

a stagger angle (Fig. 1) 

c fraction defined in equation (9) 

angular deformation (Fig. ll 
4> angular coordinate in rotating frame of reie ?nce 
0 for blade 1 

q? phase angle, equation (2) 

phase angle, equations (7) and (fl) 

CO circular frequency 


The use of stationary displacement sensors and 
high-response aerodynamic instrumentation for the tur- 
bomochlnery flutter investigation was recently doeu- 
montod [1,2]. The aerodynamic instrunrsntation in- 
cluded a blade-wake velocity transducer and case- 
mounted static ptcaaure transducora, while the dia- 
placemont inatrumontation conaiatod of two fibre- 
optics probes that sensed the light reflected from 
the rotor blade tipa. 

The advantages of stationary sensors over 
rotor-mounted senaora were domonatrated in these 
references particularly well due to the mistunod 
character of flutter. While a stationary sensor 
scans all the blades at a given position, a rotor 
mounted aonaor performs a point-moasuroment only for 
one blade. Largo variations in blade to blade 
vibratory amplitude a:id phaac wore reported from 
these measurements. 

Another adva:itago of these sensors is derived 
fium the fact that the phase differe::cc between 
blades In the rotating frame of rcfcre:Kc translates 
into a frequency shift in the stationary rofei-enoo 
frame. Thus, the prosonoo of several aerodynamic 
modes corresponding to different intorblade phase 
angles is evideneod in a spectrum of a ststionary 
sensor no the presence of several frequencies. Ex- 
pressed in engine orders, these frequencies are eas- 
ily identifiable because of the common non- 
integrality. 

While the pressure and the displacement sensors 
are similar in regard to these properties, it may be 
noted that the pressure sensors yield a continuous 
signal whereas the displacoiirant measuronKsne is essen- 
tially discrete. Because of the continuity property, 
the spectral analysis of the pressure data is more 
difficult. The method presented in [2] is only ade- 
quate in the region whore spatial unsteady pressure 
gradients sre small, typically at midchord position. 

At the leading edge, however, largo spatial variation 


in amplitudo and phase occurs within the blade pas- 
sage (3J, implying a wide frequency bandwidth of the 
stationary preoaure signal. Transformation from a 
stationary into a rotating frame of reference ouch as 
described in |2] boeomoo then impractical. 

This paper describes the alternate sampling and 
analyses techniques designed to eliminate the short- 
coming of the previous procedure. The pressure in 
each blade passage is sampled at the same position 
relative to the equilibrium blade position, as in (2) 
for the discrete blade displacement signals. This 
avoids the superposition involved in the transforma- 
tion into a rotating frame of rcforonce, but intro- 
duces another difficulty associated with the fact 
that blades are vibrating and not still. 

To deal with this problem two schemes are pre- 
sented which result in blade-p .sourc-difference 
spectra. One scheme relies o>< >ne differential- 
pressure sampling in the stoop’orossure-gradient re- 
gion of the blade pasoage. and the other relies on 
the sampling in the two adjacent regions where the 
gradients arc very small. 

The paper also includes the optical displacement 
data at two chordwise poaltiona. These data ware 
sufficient to define the motion at the blade leading 
edge. Phasing of the blade-pressuro-dif forenco spec- 
tra and the displaeenxjnt spectra at the leading edge 
indicated to which extent each nodal diameter compo- 
nent is contributing to the instability. 

In addition to the above blade-p.asoage (or blade- 
order) soinpling, the paper includeo a discussion of 
engine-order sampling in conjunction with its appli- 
cability for the on-line flutter detection. 

The detailed documentation of the test condi- 
tions, the flutter boundaries, and the flow condi- 
tions in the relative frame of reference is available 
in [4]. The test point discussed in this report is 
125. 

EXPERIMENTAL MEASUREMENTS 

The location of the instrumentation ports for the 
two types of measureroenta analyzed in this paper is 
illustrated in Fig. 1. The blade position in this 
figure includes incremental changes caused by load 
and speed which were estimated from manufacturer's 
data. Static pressures at other chordwise positions 
were also available! however, the leading edge trans- 
ducer illustrates best the sampling and analyses 
techniques described in this paper. 

The displacement data were recorded on a multi- 
channel direct-record type tape recorder at 304.8 
cm/s (120 in/s). Another FM tape was then generated 
by replaying the original tape at 7.62 cm/s (30 in/s) 
and recording the data at 304.8 cm/s (120 in/s). The 
frequency bandwidths were 2 MHz on the direct-record 
and 80 KHz on the FM. The direct-record method was 
used for the original data recording in order to 
avoid the. frequency limitation encountered on the FM 
recorder. The pressure data, however, were recorded 
directly on the FM recorder. 

Subsequent data processing was performed on a 
software-baaed Fast Fourier Transform (FFT) analyzer 
equipped with dual magnetic disc drives. The mag- 
netic discs were used fer storage and subsequent re- 
processing of digital data. This capability essen- 
tially eliminated the need for the large computer and 
the central data processing facility which were used 
in the previous [2] di ■',>lacement-dada analysis. For 
the pressure data, Che mass data storage capability 
enabled the impleme..r. ation of special sampling tech- 
niques. The ef.fective digitizing rate was 2.09 MHz 
for the displacement data and 1.31 MHz for the pres- 
sure data. These rates were achieved by greatly re- 
ducing the speed of the FM tapes on playback. 


ANALVSIS OF THE DISPLACEMENT DATA 

As can DO seen in Fig. i, instantaneous blade-tip 
position measurcmenco were made at two chordwise 
positions. I,i addition to those, _ two reference 
pulses were recorded, the IE end the 38E pulse, which 
corresponded to the blSdo pessing frc'jiency. These 
reference signals were generated by two magnetic sen- 
sors located opposite the appropriate triggers which 
wore mounted on the rotor shaft. 

In [2] ail four channels of data were digitized 
simuitancously end stored on magnetic tapes. The 
displaccmentB wore obtained by measuring the posi- 
tions of Che trailing edge pulses relative to the 
corresponding 38E reference pulses. The IE reference 
pulse provided the orientation in the relative frame 
of reference. Prior to the spectral analysis, the 
displacomonts wore ruferensod to theio respective 
averages in order to eliminate fluctuations resulting 
from blade spacing non-uniformitica. 

In this paper the data were digitized using the 
FFT analyzer which was limited Co two channels. 
Digitized were either two of the three channels com- 
prising Che two optical data channels and Che 38E 
reference pulse. The IE reference pulse was used to 
initialize the digitization and was not itself digi- 
tized. It should be noted that either reference 
pulse could be used to calculate the displacomencs. 
The 38E reference pulse has been used due to better 
signal quality. In ocher respects the analysis of 
optical data was Che same as in {2]. 

Figures 2, 3(a) and 3(b) present the overall 
spectra corresponding to the traiiing-odge sensor vs 
38E signal (relative Co which Che displacements were 
computed), the midchord sensor vs 38E signal, and the 
trailing edge sensor vs midchord sensor. A special 
procedure was followed when these data were digitized 
in order to assure that the data acquisition was ini- 
tiated with the identical IE reference pulse. This 
procedure will be discussed later in conjunction with 
the phasing of pressure and displacement spectra. 
Because the blades were sampled successively during 
each revolution, the number of samples per revolution 
was 38, and the highest detectable frequency was, 
therefore, 19E, The number of data points used in 
the spectral analysis was 2048, which approximately 
corresponded to 53.9 revolutions. The flutter fre- 
quency obtained independently from strain-gage data 
was 8.4SE. (All available strain-gage data exhibited 
the same flutter frequency.) 

In Fig. 2 the frequencies corresponding to major 
peaks are nonintegral by the same amount as the flut- 
ter frequency, i.e. , 0.45E. The nodal diameters cor- 
responding to this fraction can bo obtained by sub- 
tracting the flutter frequency in engine orders from 
each frequency associated with this fraction. The 
possible range for nodal diameters so obtained is -8 
Co 10. Beyond this range, the nodal diameters corre- 
spond to aliased frequencies which arc associated 
with fraction 0.55E. The spectral peaks correspond- 
ing Co these frequencies ore very small. 

The phase angles corresponding to major peaks 
(i.e., corresponding to dominant frequency lines as- 
sociated with major peaks) in Figs. 2 and 3 are pre- 
sented in Table 1. They were determined neglecting 
the slight misalignment of the optical ports relative 
to the blade chord (Fig. 1). Choosing the Urailing- 
edge sensor as Che base, this correction would amount 
to the addition of 11.7 degrees to the phase corre- 
sponding to Che midchord sensor. Because Che ampli- 
tudes corresponding Co this sensor are considerably 
smaller than for the trailing-edge sensor, the cor- 
rection for the phase corresponding Co Fig. 3(b), in 
the last row of Table 1 would be considerably smaller. 
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Uccaunc {!ii> ptiaae angli’a corroo Hindiiiis to a given 
ftoquoiwv 111 IflOli' I are close to each other, and the 
apoctro to l?iiso. ,! and 3 are nearly oimilor, the 
blade ri.itton con be denenbed a# a rotation about a 
ninale point, 

Dorns Fisa, 2 ami 3, the apptosimation 
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and the geometry in Fig, 1, the rotational aKio was 
located at obout of, percent ot chord trom the loading 
edge. It in aloo pooniblo to expreaa the unito in 
Fig. 2 in termn of ansular deforiaation. Aa noted in 
rhio figure one unit coriooponds to 0.117 dcgreco. 

Ad in 12], the individual blade apoctra can be 
obtained by arraitging aampled data in groupa so that 
each group includea aucccnaive saraplea from a partic- 
ular blade only. Uecauoo the oampiing for each blade 
to IE, the resulting opectra are folded into the 
range 0 to 1/2H. Figure A prooenta the flutter 
amplitudea and phase angles fur each blade as derived 
from Che trailing*edge aonsor vs 38E pulse data cor- 
responding to 64 rovolutimia. 'Ihe phase angles were 
roforonced to blade one phase angle. 

AWAbYSIS OF 'fllE PRESSURE DATA 

Engine Order Sampling 

The pToosurc data discuBBCd in this paper are 
derived from a ainglc portion of the tape which was 
digitized and atored on two magnetic diacB. Two 
channelo of data were digitized, the preoBure signal 
from the leading edge tranadueer and the IK roferonco 
oignal. There wap a total of 60 engine revolutions 
stored dioOB and on the average the number of 
points per channel per revolution was 10295. This 
provided sufficient accuracy (i.e. i/tl0295E)) to 
illustrate the different sampling methoda preaented 
in this paper. 

The main advantage ot the engine-order sampling 
ia that It enables the elimination of the ateady- 
atatc part of the proBouro aignal ao that only the 
non-integral, flutter related apeetrsl peaks remain. 

For the spectral analysis presented in this sec- 
tion, .’’ig. 5, the pressure data wore sampled 512 
times per revolution and the total number of revolu- 
tions was 16. Xhia gave a total of 8192 points which 
was the maximum number that could be processed on the 
FFT analyzer. The nonf ) uctuating part of the pres- 
sure aignal was removed by computing the average 
pressure (over 16 revolutions) for each of the 512 
points, and then subtracting this average from the 
corresponding instantaneous pressures. 

Although there were only 16 revolutions included 
in the computation of each average, the removal of 
the engine-order peaks was fairly effective. The 
most visible 38E peak was reduced in amplitude by 
about a factor of 34. The aliasing frequency for the 
data presented in Fig. 5 was 256; however, the con- 
tribution to the power spectrum beyond 128E (Fig. 

5(b)) was very ssall which insured against aliasing 
in the range from 0 to 128E (Pig. 5(a)). 

There are two non-integral frequency families in 
Fig. 5 associated with fractions 0,45 and 0.55. The 
former ia nssociated with nodal diamotors greater 
than or equal to -8 and the latter with nodal diame- 
ters less chan -8. This can be seen by expressing 
the spectral peaks in Fig. 5 as [2], 


Aj^cos[u)^(N^ + U)j/aij.)t + qi^] 


( 2 ) 


Because only positive frequencies can be detected by 
spectral analysis, for < - cug/aip, the above 
expression takes the form. 


It IB Boon now thal, the flutter f/equcncy expreoBed 
in engine orders ie being oubtracted from the pOBi- 
tive whole owmbor (-Hi) yielding the fraction 
0,55, It is notefl that the phauo angle is reflected 
about thw horizontal axia. 

It can bo shown that if all froquencica corre- 
Bponding to the aigiK’ficant peako in Fig. 5 were 
folded about the .frequency multiples of 19B, as many 
tifflca as necesBary, oo that they fall in the range of 
0 to 198, only the frequencies which arc aiao present 
in Fig. 2 would remain. It preasuro tranaducoro wore 
mounted on each blade, the nodaf diamotors corre- 
sponding to those froquencica would be the only ones 
detected, because in this case the measurement is 
spatially discrete. It follows therefore, tliat there 
are no now nodal diameters, or modes resulting from 
the spectral ana lysis of the pressure data. The ap- 
parent nodal diameters seen by the stationary proa- 
sure tranaducoro are the result of a complex pressure 
distribution in the blade channels, the description 
of which requires a large number of spatial hartnonics. 

Although the spectral results such as in Fig. 5 
are illustrative, for the purpose of quantitative 
analysia they arc not very convenient. In order to 
obtain the nodal diaroetora, it is necoBsary to fold 
froquoneieo in Fig, 5, and in order lo tranotorsi the 
reaulta into a rotating frame of refewneo it io nee- 
oBBory to perform superposition t2j. Both require 
the knowledge of phase angles for all significant 
peaks in this figure. In view of their large number 
and a wide frequency spread, this io difficult to 
aecompllBh. The wide frequency spread also implies a 
limited accuracy duo to the memory limitation of the 
spectral analyzer. 

However, the spectral results such as in Fig. 5 
can ho used for the purpose of flutter detection. In 
this ease, the engine-order aampling could be per- 
formed directly utilizing the device termed the angle 
clock (or angle encoder) [1] aa the extonial sam- 
pler. The angle clock is caaontially a frequency 
multiplier which is capable of generating an arbi- 
trary number of pulses per revolution from an input 
IE pulse. The accuracy of this device, which io now 
commercially available is comparable to that achieved 
presently through the use of tlie analog tape and 
peripheral digital storage dovicoa. The removal of 
the non-fluctuating part of the proBsurc signal es- 
sentinlly involves tne computation of the average 
which does not entail any appreciable time delay. 

An additional advantage of the engine-order sam- 
pling is that for the undorsamplod data the folding 
occurs oithor about an integral frequency or an inte- 
gral frequency plus 1/2E, depending on whether the 
number of points per revolution is oven qr odd. 
Therefore, by solocting an even number of points per 
revolution there would be only two fractions associ- 
ated with a flutter frequency even for the aliased 
sampling. For the purpose of flutter detection, a 
more important consideration than the aliafling is 
that the total number of revolutions included in the 
sample be large, wltich facilitates the officiant re- 
moval of engine order peaks. Ab will be Bhown later, 
it is also desirable that the number of samples per 
revolution is not a multiple of the number of blsdes 
on the rotor. 

Blade-Passage Sampling 

In this section an alternate sampling and an 
analysis method for the pressure data is presented. 

The sampling is performed once per blade passage at 
the same position relative to the steady-state pres- 
sure distribution. 
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Prior to the Bnmi>ling, Etier»5foie, tlu' ot«:.«uy“ 
Btoto diotribution had to be dutorminod. Ihio was 
achieved by oeparating the data recorded on disco 
into reeorda covering the period ot one revolution, 

(IB doternuned fvom the relorcneo Us puloo data. The 
pointa vitliin each revolution wore then ordered rela- 
tive to the ocart ot the revolution and the corre- 
opondiiig pointa in each roford were averaged over a 
total number of revoluti&.a. The reauiting steady- 
state distribution conaioted of IUH95 pointa, each 
repioaentlng the average of 66 revolutlonn. Fig- 
ure b llluatratcB a typleol diatribution for a blade 
paaaacc with the right-hand gradient region proaented 
enlarged. In thia figure the liorfBontal axla in given 
in torraa of number of digitized pointa relative to an 
orbitrary origin. The baaic unit of time ia thoro- 
foro 1/U0299 K), The marked pointa in Fig. 6 aub- 
divide the gradient region into equal partn each of 
which io oeven apacoa wide. Thooo pointa repreaent 
the locationa relative to the aero croaoing where the 
proBanre wan oarapled, The zero uroaoingn themae.lvoo, 
oeing common to all pointa, wore determined in ad- 
vance and atorod on a dine. 

One important consideration fur the aampling in 
the proBoure-gradient region ia the blade motion it- 
self. In torma of unita of Fig. 6 the blade dis- 
placement wao about 5. Motion of this extent intro- 
ducea a fluctuating preaaure in addition to the fluc- 
tuationa which would be recorded if the blades wore 
not vibrating. It in poooible to demonstrate thia 
effect indirectly by showing that the vibratory mo- 
tion - as derived from the preaaure data in the gra- 
dient region - dooa not compare well with the vibra- 
tory motion derived from the optical diaplaceraent 
data. 

The blade vibration ia reflected in the preanuro 
data aa a pcrioi ie displacement of the gradient por- 
tion of the proiaure diatribution relative to its 
equilibrium or otoady-atate diatribution. The prea- 
sure data can, therefore, be proceaaed quite simi- 
larly to the optical pulno data, the nuin difference 
being that the inatantancoua blade poaition ia ob- 
tained relative to the IE reference pulse rather than 
relative to a corresponding 38E pulse. The spectral 
reaultn obtained in ouch a way are preaenced in Fig. 
7(a). If there were no interference from the pres- 
sure f luctuationB, these results would bo aimilar to 
the optical displacement apectrum in Fig. 2. The 
fact that they are different indicates that there ia 
interference from ptesauro f luetuationa. Super- 
position of two signals for an i'th nodal diameter 
can be expresaed as 

A^cos j(to^ + 6)j.Nj^) t + 

+ ^l(^) ^ [H 

where the superscripts D and F denote respec- 
tively displacements and pressures. 

Thia unconventional analysia of pressure data was 
specifically performed to illustrate the displacement 
signal contamination due to pressure fluctuations. 

The sampling was performed at a constant pressure 
level corresponding to point 7 in Fig. 6, More con- 
ventional results in Fig. 7(b), however, were ob- 
tained by reversing this process, i.e., by sampling 
preasuros at the position defined by this point. It 
is seen that results in these two figures are quite 
similar. Both represent the summation of two sig- 
nals, one derived from the pressure fluctuation and 
the other from the blade motion as shown in equation 
(4). Corresponding expression for Fig. 7(b) would 


differ only by a multiplicative taetur dP/dt. It 
would illustrate contamination of preauurc fluctua- 
tions duo to blade vibratory motion. 

(iradient Sar tiling 

‘“^OnF’aoTution to this problem ia to form a proa- 
sure difference in the region where the ntnady-atatc 
prtasure gradient ia approximately conotant. liecauoc 
the blade diaplooemont io the anmo on both aidoo of 
the blade, and the gradient ia approximately con- 
stant, this contribution to the preaaure difference 
from the blade motion ia nearly eancclled. 

To illuacrato thia, prCBouroa wore aamplcd at 
indicated pointa in Fig. 6 and then the following 
diffcrcnceo wore formed, Pj - Pq, P^q - Py, 

Pl 4 ” Pq end P 2 f - P-y. The opeetrii corro- 
npotiding to thsou difforencon are given in Fig. g and 
phase angle, are prcaonted in Table 2. The number of 
aamplcd data points used for apectral analysis wao 
2048. Becauao the presnorcs in each difference above 
wore not nampled at the name time, a more exact pro- 
cedure would be to obtain tranaforraa for each point 
individually, apply the appropriavc correction to the 
phase angle of one point and then perform the nub- 
traction in the frequency domain. However, the maxi- 
mum correction for the phaao angle tor the reaulta 
preoonted in Fig. 8 would be 8 degrooo. Therefore, 
the reaulta obtained uoing the latter procedure would 
not differ appreciably from the above. 

It can be seen that the prcaauro-differeneu spec- 
tra in Fig. 8 are all qualitatively similar and that 
they are oimilar to tlie displacement apectrum in 
Fig, 2. Comparison of the aroplitodes of correa-ind- 
ing puaka in Fig, 8 aiao indicatea that in the ap- 
proximately constant-gradient region between pointa 0 
and 14 IFig, 6) the amplitudes arc approximately pro- 
portional to the number of intervals between the pair 
of points included in the difference. If thia intei- 
val is maintained constant the exact poaition of the 
pair of points does not seem to bo important. Com- 
parison of phase anglea at a given frequency in Table 
2 indicatea that they are noariy constant. 

These observations can be explained qualitatively 
based on physical arguments. Because it takea a 
finite time to traverse the static preaaure port and 
because of the finite blade thickneas in the direc- 
tion of the wheel speed, it is expected that the 
pressure difference noroae the blade will be smeared 
over approximately 10 poinls. This corresponds ap- 
proximately to half the width of the ateep gradient 
region in Fig. 6. It appears, therefore, that the 
smearing of the preeaure difference is in part duo to 
dimensional effects. 

Considering first the smearing due to static port 
dimensions only, one can treat the blade aa a dio- 
continuity as shown in Fig. 9, The sketch on the 
left (‘orrosponds to the left-point gradient sampling 
and on the right to the right-point gradient sampling 
(Fig. 6). The inatantanooua blade position ia de- 
fined in terms of the area fraction f, where f = 0 
corresponds to the position of the blade when it ia 
aligned with the left edge of the port. Subscripts 
L and R denote left and right positions, and P 
and S denote the blade pressure and suction sides. 
The pressure difference between tlie right and left 
pressures is then 

Pr - “ Vf - '^r)’’s ■ fVp 

or, 

^R * ’^L “ ' "^P^ 
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tiulCP aU i'U'uauK'B au’ vcvti'ioi poBaospinfl ao” 
pljcuijp onJ [ilmai’i ami i jti 's'ni, it tuiluwa that 
tliO ritt'onure liitlcnnuT hi'twca the twe pointo in the 
conBtant“gradinn rciUim haa tlU' oaat> pliaoc aa the 
auk tion =f>ioaBiitr aulc ilittincncp. It can alBik bo 
aocn that ttio r'ik'P>nUcnaUty i»t oiriUitudo am) tlio 
diotaucp betworn the twn iiouita ahouli! lioli) oppioxi- 
oately in the cential tegion whoro Itactivin t la 
apiiroKiaiatolv pu'pk'ii ik'nal tik the blade poaitien. 

tiainp, two appri'Mr-iati.mai an eKplanatien Ckinld 
alae be ntteied tor tlie eaae et Unite blade thick* 
noan an*) nrslipibU' pteanuie poit dineuoik)na< The 
Iirat 10 the quani ateadyatnto aaauniption, baaed on 
the tact that the blade thicknoaa ie araall corapaicd 
to the wavolenath conuBpoudiii>i to one tlutter cyclOi 
and the aecond la the linoarily of picaflure variation 
in the tip»cloarance tegion over the blade thicknoaa. 

2e/>'_^;imlient.J)^ 

Uccavtae the urteitetence with preaauto lluctua* 
tiona doe to blade notion lo proportional to the 
ateady-atatc pieaaure pradienti an alternate way ol 
eliminating thia interterunee is to aaraplo in the 
legion where preaauro Biadicntn arc amall, or ideal- 
ly, sero. For example, Fig. ItKa) preaentn the opec" 
trum tor point - Id (Ftp,, b) whleh ia atlll tnlrly 
elone to the atcep grauient region of clic proaaure 
diatribution. It can be aecn that thia apeetrum ia 
qualitatively much cloaer to the diaplacemont opoc- 
tvum in Fip,, J than the one in Fig. Kb), which waa 
obtained by aaiiipling in the ateep gradient portion of 
the blade pannage. 

In order to obtain the preonurc dillerencc acroaa 
the blade, one aloo baa to aeleet a point in the zero 
alopc region on the auction aide Of the blade. Uc" 
eauoe ot comudciable curvature in the ateady-atatc 
proBauie Uiotribution in tlua region and the blade to 
blade diifeteucea, tlii« ia more diltieult. lianed on 
roaulta obtained oo inr, one would expect that the 
opeccrum lor thio point oliould also be qualitatively 
aimilar to the diiiplacement Bpectrum in Fig. 2. It 
wao found that thio condition ia ronaonahly oatiofiod 
for point 2ii IFig. b). The apeetrum lor thiu point 
ifl preaentod in Fig. llHb), the preaoure difference 
apeetrum corro.nponding to pointa 28 and -14 in Fig. 
10(c), and ita phaue in Table 2. TTie apeetral ampli- 
tude and phaae for the preaoure difference IMg - 
I’-iq if! ueen to be cloae to the reopeciive vafuea 
obtained lor I’ji " l’- 7 > 

Although it ban not been atated explicitly tliua 
far, the underlying nooumptiun tor obtaining a valid 
preaoure difference ucroao the blade io that the, 
preoBure variation over an interval correopoiiding to 
the diotanee between the two pointa oolccted tor com- 
putation ol the dilforence ahould bo small, except 
when erooaing the blade. One way to chock this aa- 
sumption is to fom a preasure diltcrunce in the ap- 
proximately flat regioo to the left of the zero 
croBsing in Fig. b, e.g. , F-? - It io not 

likely tlian those two poinds will bracket the bJadc 
and therefore one should expect amplitudes which are 
much amallor than a difference involvitig the same 
number of points in the gradient region. This is 
subBcantiatad by comparing results in Figa. lOtd), 

8(a) and 8(b) corresponding roopeetively to - 
l’-14, 1’7 " I’o l’l4 “ )’? difference. 

Smear ing Ef fects Du e to bla de Mot ion 

iiT a previous ieetion tiie pressure-difference 
smearing due to the finite dimensions of the trans- 
ducer and the blade was examined. However, even in 
the absence of dimensional effects, in tlie gradient 
region some smearing can also occur duo to the blade 
motion. Because the two sampling positions for the 
pressure-dif teronce spectra were selected bnand on 


the average preasuio dintribution, and becaudo ol 
blade vibratory »tion, it lo ponoiblo to have both 
prcBSute aamploB fotning a dtflerenco on the aamo 
side of the blade at least soractimoo. Thro would 
tend to ted lie 0 the magnitude of the preanuro differ- 
ence. 

To srmnlate thia effect, two data sotB are uti- 
lized aaoociatcd rcapeetivoly with the opectra in 
Figa. 10(d) and 8(d). The first act, correoponding 
to Fig. lOCd) io roprcaontativo of tlie prenouie dif- 
ference when both preosuoo arc on the oamo aide of 
the blade and the accoK',^, corrooponding to Fig. (Hd), 
is representative of preaauro ditleronee data when 
the two piesauroB fonning the difference are on the 
oppoBito oidea of the blade. A third data set was 
then created by random sampling of the firat data ooc 
followed by oubstitution of the sampled valucn into 
corresponding point® the aecond set. Only pointa 
which huve not aeon aelected pruviouBly were sub- 
otituted in cho aecond act. In the particular case 
considered, the third data act conaiated of 1024 
points from the first and 1024 points from the second 
sot. Spectral amplitiidoa fur thcae data are pre- 
sented in Pig. U and phase angles in Table 1. It 
can bo seen that the apoctrai peaks are approximately 
proportional to the number of pointn from the data 
sot eorreoponding to Fig. 8(d), and that phaae angleo 
agree fairly well with other corroiipmiding phaoo 
angleo in Table I, Conoidering the magnitude of 
poako in Fig. U one can aooociate thio apeetrum with 
the one in Fig. 8(c) which corteoponda to the pres- 
sure difference I'lq - Pg. Qualitatively, one of 
the principal difforoncea between thoao two figuwo 
ia the Uigbnr noise level in Pig. U. However, thio 
may bo due to the fact that the pressure difference 
in Fig. U originates from two pointa which are apart 
twice the distance involved in Fig. 8(c). The larger 
distance between the points in Fig. 11 implies that 
the noise between those two points is less correlated 
and that, therefore, during aubtractiou the cancella- 
tion of noise is loss pronounced. 

CORRHbATtON OF PRESSURE AND DISPbACEHENTS 

111 order to be able to correlate the measurements 
recorded on diflerent tapes it was essential to be 
able to initialize digitization at exactly the same 
time. Thia was achieved by using a time code reader 
which produced a pulse at the instant the continuous- 
ly rend time from the tape matched the preselected 
time dialed on the time code reader. The data ac- 
quisition was initiated at the first IE reference 
pulse following this instant. The coincidence condi- 
tion wan detected under software control and the IE 
reference pulse was used as the external trigger. 

The suggestion Cor the use of the time code render 
and its implementation in thia procedure was made by 
H, It, Jones of NASA bowis. 

Because the time code for the diaplaaement data 
was not available on the original d’ '•'c-t-rccord tape, 
it was not possible to corielato directly the pres- 
sure and displucement data. (An arbitrary referenced 
time code, however, was recorded on the FH copy of 
the dlsplnoeiiK’nt data tape enabling repented replay 
of this tape at the identical IE reference pulae). 

In order to correlate tlie displacement and the pres- 
sure data, it was necessary, tborefore, to use the 
strain-gage data for which the time code was avail- 
able. Thia was accomplished using a procedure which 
is best illustrated by reviewing the resulting ex- 
pressions for the correlated displacement and prea- 
sure signals. For an arbitrary i-th nodal diameter 
those expressions are 
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P,D » Aj'®oiip j[wt - (« - + 5^'®] (6) 

1 .1 

vhcrc P »."< D denote roopectivcly prcoourc and 
dioplaeeiaont. t\t.‘ etencc time t “ 0 corrosponda to 
the isiotant a particular IIS roierenco puloo roacheo a 
ptedetorrained level and, ao indicated, tlie reference 
angle corteopoitda to blade 1 angular pooitiun. Iho 
phaoc angloa ^'i for the dioplaecments and proa- 
o.ireo are given rcopectively ao 



In these expressions tpj,'o are the i-th nodal diam- 
eter phase angles moasured relative to the first 
point in the input data set for spectral analysis. 

For the displacement data, this point corresponded to 
the first blade that appeared under the sensor fol- 
lQ«ing the IE trigger and for the pressure data it 
corresponded to the second blade after the IE trig- 
ger. (The latter choice was necessary because the 
first blade under the sensor was slightly past the 
selected triggering level at the time of IE trig- 
ger.) The second terms in those expressions account 
for the position of blade one relative to the above 
reference blades for <Pi's, and Che third terms 
account for Che time delay between the occurrence of 
the IE pulse and the instance appropriate reference 
blades are under the respective sensors. The fourth 
term in the expression for $5, <p| is Che phase 
angle for the blade-one above-shround-traiiing-edge 
gage (ASTE gage [4]) relative to the IE trigger, and 
the last term accounts for the arbitrary time refer- 
ence for the displacement data. This terra represents 
the phase of blade one at the ( imo of IE trigger. It 
was actually evaluated from th< aliased 1/2E spectrum 
associated with this blade, and not from the given 
expression. 

The strain-gage phase vf, was determined by 
interpolation between the frequency lines, similar as 
in [2]. However, a simpler and more accurate approx- 
imation was used for fraction e which defines the 
true frequency relative to a neighboring frequency 
line. This approximation is due to D, Braun of NASA 
Lewis. If cUf is the true frequency and the two 
most significant amplitudes are Aj and A2 
corresponding to frequencies u)]^ and 0)2, and 
(JUI < CDjf < CP2, then 

(‘“f ■ ^2 

^1 " 2tt " + A^ 

The phase angle corresponding Co cuf is equal to 
the phase angle at CU]^ minus eiTr(n - l)/n ~ eiTT. 

In practice, the accuracy with which eq can be 
determined la usually limited by the signal to 
noise ratio, An estlraaCo on the accuracy can 
therefore he made using this ratio, Another esti- 
mate on the accuracy may be obtained by varying 
parameters in the spectral analysis. Sampling in 


the range of about 8 to about polnto per cycle, 

0 spread in tlio phase angle of aheut 10 dc» 
greco was encountered. It is therefore eotlraated 
that this phase angle lo accurate within -r^ depreea. 

Because of the avaxlabiixcy of the liingletcn 
FFT-algorichm [5], the spectral analysis ot the dio- 
placoment and prosouro data discussed in this oei.tiuu 
was performed on a largo computer. This algutithrj 
does not require that the number ot pointa for opcc- 
tral analysis be a power of two. It is poaoihle, 
therefore, to select the aite of the input data block 
for the spectral analysis consisting of an integral 
number of revolutions and having approximately an 
integral number of cycles. This results in a unifonr. 
nonintegraiity of flutter peaks and minimizes the 
fraction cj^ (computed relative to the frequency 
line corresponding to the dominant peak) , 

For the displacement data, the number ot data 
points was 2470, corresponding to 65 rovolutionu and 
containing 549.1 flutter cycles (aosoeioted with 
Ni “ 0). For the presaute data, the number of data 
points corresponded to 63 revolutions, and the number 
of flutter cycles was 523,9, Different choice for 
the number of revolutions for the pressure data was 
governed by a alight ahift in the flutter frequency 
expressed in engine oidora. 

The phase angles corresponding to different peaks 
in the respective spectra were obtained by using an 
effective fraction e. For the displacement phase 
angles, it was obtained by averaging the c's corre- 
sponding to six most significant peaks, and for the 
pressure phase angles, it was obtained by averaging 
values corresponding to two most significant peaks. 
The phase angias were estimated to be accurate within 
about ±10 dogroea, 

The phase difference between prassure and dis- 
placement vectors is presented in Table 3. lieeauoo 
pressure loads the displacement, all modes in this 
table are unstable. 

It is noted that precise determination of the 
moment coefficient and the work per cycle vrauld re- 
quire a repeated application of the methods outlined 
in this paper at the remaining chordwiso locationa. 
However, because the unsteady force on the blade dur- 
ing Stoll flutter is concentrated near the leading 
edge [3,6], the conclusion regarding the stability ot 
different modes reached in this paper is not expected 
to change by this extension, 

DISCUSSION OF RESULTS 

In reviewing the different method of anolyaio of 
pressure data presented in this paper, it may hu 
noted that basic to all is the engine-order sampling 
followed by the removal of the steady-state part of 
the pressure signal at each sampled point. If pres- 
sure signals are processed mainly for flutter detec- 
tion, as already mentioned, it is desirable that the 
number of samples per revolution be not an integral 
multiple of the number of blades. This assures that 
pressures at different relative positions in the 
blade passage are sampled ao that the possibility of 
sampling only in a region where the unsteady pressure 
signal is small, or perhaps not even detectable, is 
eliminated. For example, sampling once per blade 
passage at a poind located at 110 intei'valo relative 
to the zero crossing (Fig. 6) results in a spectrum 
with 11.45 peak reduced by a factor of 4.6 and the 
5.45 peak by 0 factor of 2.3 relative to the corre- 
sponding peaks in Fig. 10(b). 

Blade-order sampling is therefore better suited 
for post-run data analysis when the precise blade- 
passage position of the sampled point can be con- 
trolled and the possible interference duo to the 
steady-state pressure gradient taken into account. 
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Two ncthojB n»o available to accojaplioh this, tlio 
gradioiiC pK>ooaie«diffetenco oatspling and the s«ro-» 
gradient oaapiing. The toiDor method can only pro- 
vide preoaiiu'-dilEorenee apectra while the latter 
a loo violaa the prooauto opeetra cloao to the blade 
sutlacos, Ingv^noral, tlio applieability of c./ch 
Bothod will depond on the eliaractor of the steady- 
atate dietribution. In the particular eaoo, the 
gradient method appearo to be easier to apply. Thio 
Method yielda correct phase reoulta, and the arapli- 
f idea are attenuated approximately an the ratio of 
the distance between the two pointo to the distance 
correoponding to the full extent ot the gradient re- 
gion, 

It ia noted that in the computation of the un- 
steady work per cycle, the intortcrcnec due to the 
blade motion is of no oonaoquence since it in either 
in or out ot pbano with the dioplacemonc so that it 
does not. contribute to the work integrated over a 
flutter cycle. 

The analysis ot preaoure data has been give" thua 
far in termn of the over-all spectra. An anal > mn 
procedure would be to use the individual blade spec- 
tra, sitallnrly to Fig. 3 for the displaccmenta, Fig- 
ure 12, for example, illustraton auch rcaults for the 
presoure difference - P-;, Such reoulta 
were also obtained for all pointa in Fig. b and oil 
the differential prcaauroa aaoociaced with theae 
pointa. However, the over-all apectra, which exhibit 
aystem chsractorintico, wore found to be more con- 
venient to uoe for the anaiyoia. Reoulta ouch ao in 
Fig, 12 load to a diocuaoion of the individual blade 
charactoriotica whicli appeared to bo loan consiatont 
than the ayutem aa a whole, 

CONUhUDINi; REMARKS 

An analyoio procedure baoed on engine-order oam- 
pling waa prooontud which effectively cliroinatoo the 
engine harroonico from the overall opeetra corro- 
oponding to a otationary preooure tranoducer. It woo 
ohown how thio procedure can he implemented for the 
on-line flutter detection, 

Quantitative opectral annlyooo of proaouro data 
were perfotroed baaed on blade order oarapling. Opti- 
cal dloplacement opcctro wore uaed to demonstrate the 
intorferonco of blade motion with the preooure signal 
in the atoop gradient portion of the blade paoooge. 
Two metbode were outlined which yield tha blade- 
proaouce-differeneo opcctra removing the contribution 
due to the blade motion, 

Ph.oaing of the blade differential preaoure and 
the diaplaecraent opeetra at the leading edge indi- 
cated that preooure loads tin? dioplacoment for the 
six moot sir- ficont nodal-diomoter inodes, 
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TABLE 1. - PHASE ANGLES, DISPUCEMENT SPECTRA 


Figure 

Frequency (B) 


3. AS 

4. AS 

5.45 

8.45 

11.45 

12.45 

13. t5 

2 

-65 

-61 

-57 

48 

-68 

-68 

-65 

3(a) 

-56 

-SA 

-49 

77 

-83 

-79 

-86 

3(b) 

-66 

-64 

-59 

46 

-63 

-65 

-56 


TABLE 2. - PHASE ANGLES, PRESSURE-DIFFERF.NGF SPECTRA 


Figure 

Frequency (E) 


3.45 

4.45 

5,45 

8.45 

11.45 

12.45 

13,45 

8(e) 

-144 

133 

43 

-132 

-161 

98 

28 

8(b) 

-144 

129 

59 

-134 

-161 

101 

38 

8(e) 

-144 

131 

51 

-133 

-161 

99 

33 

8(d) 

-149 

127 

47 

■ 132 

-163 

90 

36 

10(c) 

-142 

131 

48 

-136 

-171 

90 

27 

U 

-133 

128 

48 

-136 

-164 

95 

35 




TABLE 3. - DIFFERENTIAL-PRESSURE PHASE REUTIVE 
TO DISPUCEMENT PHASE 


AP 

Nodal diameter 

-5 

-4 

-3 

0 

3 

4 

5 

^4->’0 

7 

11 

14 

25 

33 

24 

51 
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i ,^AXIAL POSITION OF THE 

LEADING-EDGE PRESSURE PORT 



I 

Figure 1. - Location of measurements ports. 
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Figure 6. - Stra<^-stale pressure distibulion; 1 unit 
mmsB. 















